Measurement of e+e D^*^^ D^*'^ cross sections near threshold using initial-state 

radiation 

G. Pakhlova/^ I. Adachi/o H. Aihara.^o K. Arinstein,!- T. Aushev,24. 17 t. Aziz^s A. M. Bakich,^^ 
V. Balagura,!^ E. Barberio,^^ A. Bay,^^ K. Belous/^ V. Bhardwaj,'*° B. Bhuyan,!^ A. Bondar/'^^ 
A. Bozek,34 M. Bracko,^^.!^ t. E. Browder,^ A. Chen,3i P. Chen,33 B. G. Cheon,^ R. Chistov/^ I.-S. Cho,'^^ 
K. Cho,2i K.-S. Choi,'55 S.-K. Choi/ Y. Choi,^^ j_ Dalseno,27.46 m. Danilov,!^ Z. Dolezal,^ A. Drutskoy,^ 
S. Eidelinan,i'38 D. Epifanov/^^s Gabyshcv,!' A. Garmash/'^s B. Golob,^^. is H. Ha,22 J. Haba/^ 
K. Hayasaka,29 H. Hayashii,30 Y. Horii,49 Y. Hoshi,*^ W.-S. Hou,^^ H. J. Hyun,23 T. lijima,^^ K. Inanii,^^ 
R. Itoh/° M. Iwabuchi,55 Y. Iwasaki/o N. J. Joshi,45 T. Julius,^^ J. H. Kang,^^ P. Kapusta,^^ H. Kawai,^ 
T. Kawasaki,36 C. Kicsling,^^ H. J. Kim,23 H. O. Kini,23 M. J. Kim,23 S. K. Kim,^^ Y. J. Kini,^ K. Kinoshita," 
B. R. Ko,22 S. Korpar,26^i8 P. Krizan.^^.is T. Kumita,^i A. Kuzmiii,i'38 Y.-J. Kwoii,^^ S.-H. Kycong,'^^ 
J. S. Lange,5 M. J. Lee,''^ S.-H. Lee,^^ C. Liu,^! Y. Lm,33 D. Livcntsev/^ R. Louvot,^^ A. Matyja,34 S. McOnie ^4 
K. Miyabayashi,30 H. Miyata,36 Y. Miyazaki,^^ R. Mizuk/^ G. B. Mohanty,^^ T. Mori,^^ Y. Nagasaka," 
E. Nakano,39 M. Nakao/° H. Nakazawa,3i S. Nishida/o K. Nishiniura,^ O. Nitoh,52 T. Nozaki/" S. Ogawa,^^ 
T. Ohshima,29 S. Okuno,i^ S. L. 01sen,42.9 p. Pakhlov,!^ H. Palka,34 C. W. Park,43 H. Park,23 H. K. Park,23 
R. Pcstotnik,i8 M. Petric/* L. E. Piilonen,53 A. Poluektov,i'38 m. Rohrken,^^ S. Ryu,^^ H. Sahoo,^ K. Sakai,io 
Y. Sakai/o O. Schneider,^* C. Schwanda,!^ K. Senyo,^^ M. E. Sevior,^^ M. Shapkin,!^ V. Shebalin/^ 38 q p ^^^^^9 
J.-G. Shiu,33 B. Shwartz,i'38 F. Simon,27.46 Smerkol,!^ Y.-S. Sohn,^^ A. Sokolov,!^ E. Solovieva,!^ S. Stanic,37 
M. Staric,!® T. Sumiyoshi,^! Y. Teramoto,39 I. Tikhomirov,!^ K. Trabelsi,i° S. Uehara,i° T. Uglov,i^ Y. Unno,^ 
S. Uno,i° G. Vamer,9 K. E. Varvell,^* A. Vinokurova,i'38 a. Vossen,!^ C. H. Wang,32 ^ Wang,33 
P. Wang," M. Watanabe,36 Y. Watanabe,!^ R. Wedd,^^ E. Won,22 Y. Yamashita,35 C. Z. Yiian/^ 
Z. P. Zhang,4i V. Zhilich,i'38 P. Zhou,^^ V. Zhulanov,i'38 T. Zivko/^ A. Zupanc,^^ and O. Zyukovai'38 

(The Belle Collaboration) 

' Budker Institute of Nuclear Physics, Novosibirsk 
^Faculty of Mathematics and Physics, Charles University, Prague 
Chiha University, Chiha 
^ University of Cincinnati, Cincinnati, Ohio 45221 
^ Justus-Liebig-Universitdt Gieflen, Ciefien 
^The Graduate University for Advanced Studies, Hayama 
^ Cyeongsang National University, Chinju 
^Hanyang University, Seoul 
^University of Hawaii, Honolulu, Hawaii 96822 
^"High Energy Accelerator Research Organization (KEK), Tsukuba 
Hiroshima Institute of Technology, Hiroshima 
''^University of Illinois at Urbana- Champaign, Urbana, Illinois 61801 
'^^ Indian Institute of Technology Guwahati, Cuwahati 
''^Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 
'^Institute of High Energy Physics, Vienna 
Institute of High Energy Physics, Protvino 
'''institute for Theoretical and Experimental Physics, Moscow 
^*J. Stefan Institute, Ljubljana 
''^ Kanagawa University, Yokohama 
^''institut fiir Experimentelle Kemphysik, Karlsruher Institut fiir Technologic, Karlsruhe 
Korea Institute of Science and Technology Information, Daejeon 
'^'^ Korea University, Seoul 
^^Kyungpook National University, Taegu 
Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne 
'^^ Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 
University of Maribor, Maribor 
Max- Planck- Institut fiir Physik, Miinchen 
'"^^ University of Melbourne, School of Physics, Victoria 3010 
^''Nagoya University, Nagoya 
^"Nara Women's University, Nara 
^'National Central University, Chung-li 
"''^National United University, Miao Li 



2 



Department of Physics, National Taiwan University, Taipei 
Niewodniczanski Institute of Nuclear Physics, Krakow 
'^^Nippon Dental University, Niigata 
''^Niigata University, Niigata 
''"^ University of Nova Gorica, Nova Gonca 
^^Novosibirsk State University, Novosibirsk 
Osaka City University, Osaka 
^'^Panjab University, Ghandigarh 
'^^ University of Science and Technology of Ghma, Hefei 
^^Seoul National University, Seoul 
^'^ Sungkyunkwan University, Suwon 
School of Physics, University of Sydney, NSW 2006 
Tata Institute of Fundamental Research, Mumbai 
'^''Excellence Gluster Universe, Technische Universitdt Miinchen, Garching 
'^^Toho University, Funabashi 
'^^Tohoku Gakum University, Tagajo 
Tohoku University, Sendai 
^'^ Department of Physics, University of Tokyo, Tokyo 
Tokyo Metropolitan University, Tokyo 
^"^ Tokyo University of Agriculture and Technology, Tokyo 
^'^CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 
Wayne State University, Detroit, Michigan 48202 
Yonsei University, Seoul 

We report a measurement of exclusive e"^e~ D^*'''^ Di*''~ cross sections as a function of center- 
of-mass energy near oi*-''^ oi*-' threshold with initial-state radiation. The analysis is based on a 
data sample collected with the Belle detector with an integrated luminosity of 967 fb~^. 
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Recently a number of measurements of exclusive cross 
sections for e'^e~ annihilation into charmed hadron pairs 
above open-charm threshold were performed by the B- 
factory experiments using initial-state radiation (ISR). 
These include Belle measurements [l| of e+e" cross sec- 
tions for the DD {D = D° or D+), D+D*-, D*+D*-, 
D°D-Tr+, L>°£)*-7r+ and A+A7_final_states [|Ji and 
BaBar measurements of the DD, DD*, D*D* final 
states 0,11], which are, in general, consistent with those 
of Belle. In addition, CLEO scanned the e'^e~ center of 
mass energy range from 3.97 to 4.26 GeV and obtained 
exchisive cross sections for the DD, DD*, D*D*, DDtt 
and DD*TT final states 0]. The first measurements of 
the exclusive cross sections for e^e~ annihilation into 
charmed strange meson pairs Z?i*^^£)i*^ were performed 
by CLEO with high accuracy but with limited maximum 
energy (4.26 GeV) ^ and by BaBar ^ with a 100 MeV 
bin size. The observed cross sections were found to be 
an order of magnitude smaller than those for non-strange 
charmed meson production. 

Although the recent BES fit to the total cross section 
for hadron production in e+e^ provided new parame- 
ter values for the i/j resonances the available exclu- 
sive e+e" cross sections have not yet been qualitatively 
explained. One of the main problems is the numerous 
open charm thresholds in the region that influence the 
cross section behavior and, thus, complicate theoretical 
descriptions. 

The Y states [l^l (with masses above open charm 
threshold and quantum numbers J^^ = 1 ). which do 



not exhibit strong decays to any of the measured open 
charm final states and have remain unexplained since 
their discovery more than five years, provide additional 
motivation to pursue all possible experimental informa- 
tion about the decomposition of charmed particle pro- 
duction in the charm-threshold region. 

Here we report a measurement of exclusive e+e^ 

cross sections as a function of center-of-mass 
energy from the Ds Dl thresholds to 5.0 GeV, con- 
tinuing our studies of the exclusive open charm produc- 
tion in this mass range. The analysis is based on a study 
of events with ISR photons in a data sample collected 
with the Belle detector at the T{nS) [n = 1,...,5) res- 
onances and nearby continuum with an integrated lumi- 
nosity of 967 fb~^ at the KEKB asymmetric-energy e+e" 
collider. 

We follow the full reconstruction method that was pre- 
viously used for the measurements of the e^e~ cross sec- 
tions to DT), D°D-tt+ and D'^D*-n+ final states (2, 

13, Hi- We select e+e^ Di*^^ Di*^~^isK signal events 
in which the -Oi*'^ and ui*'' mesons are fully recon- 
structed. In general, the 7isr is not required to be 
detected and its presence in the event is inferred from 
a peak at zero in the spectrum of recoil mass squared 
against the Ds Ds system. The recoil mass squared 
is defined as: 
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where i?c.m. is the initial e+e center-of-mass (cm.) 



energy, E 



and 



are energy and 



three- momentum of the Ds Ds"' combination, respec- 
tively. To suppress backgrounds two cases arc consid- 
ered: (1) the 71SR is outside of the detector acceptance 

and the polar angle for the Ds Dl combination 
in the cm. frame is in the range | cos(0^(*)+^(.)-)| > 

0.9; (2) the 71SR is within the detector acceptance 
(I cos(0^(.)+^(.)-)| < 0.9). In the latter case, the 
71SR is required to be detected and the mass of the 
71SR, combination should to be greater than 



(^'C.Il 



0.5 GeV). To suppress backgrounds from 
e+e^ — Di*^^I?i*^~(n)(7r+7r~)7isR, (n > 0) processes, 
we exclude events that contain additional charged tracks 
that were not used in the Dl ' and the Dl ' recon- 
struction. 

All charged tracks are required to originate from the 
vicinity of the interaction point (IP); we impose the re- 
quirements \dr\ < 1cm and \dz\ < 4cm, where dr and 
dz are the impact parameters perpendicular to and along 
the beam direction with respect to the IP, respectively. 
Charged kaons are required to have a ratio of particle 
identification likelihood, Vk = I {C.K'V L-,^) > 0.6 [l3| . 
No identification requirements are applied for pion can- 
didates. 

Kg candidates arc reconstructed from tt+tt^ pairs with 
an invariant mass within 10 McV/c^ of the Kg mass. The 
distance between the two pion tracks at the Kg vertex 
must be less than 1 cm, the transverse flight distance 
from the interaction point is required to be greater than 
0.1 cm, and the angle between the Kg momentum direc- 
tion and the flight direction in the x — y plane should be 
smaller than 0.1 rad. 

Photons are reconstructed in the electromagnetic 
calorimeter as showers with energies greater than 50 MeV 
that arc not associated with charged tracks. Pairs of pho- 
tons are combined to form 7r° candidates. If the mass of 
a 77 pair lies within 15MeV/c^ of the tt" mass, the pair 
is fitted with a tt^ mass constraint and considered as a tt" 
candidate. ISR photon candidates are required to have 
energies greater than 2.5 GeV. Photon candidates used 
in 77, rj' and D*^ reconstruction arc required to have en- 
ergies greater than 100 MeV. 

77 candidates are reconstructed using tt^tt~tt'^ 
(ilOMeV/c^ mass window) and 77 (±20MeV/c2 
mass window) decay modes (~ 2.5 cr in each case), rj' 
candidates are reconstructed using •qn^-K^ (±10MeV/c^ 
mass window) and jtt^tt~ (±15MeV/c^ mass window) 
decay modes (~ 2.0 cr in each case). A mass- and 
vertex-constrained fit is applied to 77 and rj' candidates. 

D'^ candidates are reconstructed using six decay 
modes: K^gK+ , K-K+tt+, K-K+tt+tt^, K^gK~TT+TT+, 
T]TT^ and 7/'7r+. Before calculation of the candidate 
mass, a vertex fit to a common vertex is performed for 
tracks that form the Df candidate. A ±15 MeV/c^ mass 
signal window is used for all modes (~ 3cr in each case). 
To improve the momentum resolution of meson can- 



didates, the tracks from the -D+ candidate are fitted to 
a common vertex with a mass Df mass constraint. D*~^ 
candidates are reconstructed using the Df'f decay mode. 
A ±15MeV/c^ mass window is used (~ 2.5 cr). A mass- 
constrained fit is also applied to candidates. 

The Df and sidebands used for background stud- 
ies arc four times as large as the signal region and are di- 
vided into windows of the same width as that of the sig- 
nal. To avoid signal over-subtraction, the selected side- 
bands are shifted by 30MeV/c^ from the signal region. 
The Ds{D*) candidates from these sidebands are refitted 
to the central mass value of each window. 

The M^^^^-^{D+Dj) distribution for Mj^+jy- < 
5.0 GcV/c^ after all the requirements described above 
is shown in Fig. [T]a). A clear peak corresponding to 



the e"* 



Dj 71SR process is evident near zero 



recoil mass. The shoulder at positive values is mainly 
due to e+e~ — >■ I?+£'*~7isr background. We de- 
fine the signal region with an asymmetric requirement 
-0.7(GeV/c2)2 < Ml^^,,{DtD:) < OAiGeY/c^ to 
suppress event this background. The polar angle distri- 
bution of the Df Dj combinations and the mass spec- 
trum of the DfDjjisYi combinations (after subtraction 
of Ec.m ) in case (2), after M^^^^^^{Df Dj) requirement, 
are shown in Figs. [T]b,c). These distributions are in 
agreement with a MC simulation and are typical of ISR 
production. The A'/^+^- spectrum with all the require- 
ments applied is shown in Fig. [T]d). A clear peak is seen 
at threshold near the -0(4040) mass. 



L a) 









r c) 











1.5 -1 -0.5 0.5 

M(D/D,-y,,«)-E^,„,GeV 




M(D +D -), GeV/c 



FIG. 1: a) The distribution oi M^^^^n{Dt Ds) for M^+^- < 
5.0GeV/c^ after all the requirements are applied, b) The 
polar angle distribution of the D~ combinations, c) The 
mass spectrum of the DtDj^isB. combinations after subtrac- 
tion of -Bern, energy in case (2). d) The M^+^_ spectrum 
after all the requirements applied. Cross-hached histograms 
show the normalized Af^+ and A/^- sideband contributions. 

Feed down from the DfDg^ final state is shown by the open 
histograms. The signal windows are shown by vertical dashed 
lines. 

The Ml^^.^{DfD;-) distribution for Mjj+^j,- < 
5.0 GeV/c^ after all the requirements are applied is shown 
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in Fig. [2]a). A clear peak corresponding to e+e~ 
£'+D*~7iSR signal process is again evident around zero. 
We define the signal region for AI^^^^^^{Df Dl~) by a re- 
quirement ±0.7( GeV/c^)^ around zero. The polar angle 
distribution of the DfD*~ combinations and the mass 
spectrum oi Df D*^ jisn combinations (after subtraction 
of the Ec.m. energy) in case (2) after the requirement on 
M^ecoiii^t^t") is applied are shown in Figs.[2]b, c). The 
Mjy+jyt- spectrum after all the requirements are applied 
is shown in Fig. [2]d). Two clear peaks are seen at the 
V'(4160) and the 7/)(4415) masses. 
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FIG 



a) The distribution of the M^^^„n{Df D* ) for 
< 5.0GeV/c^ after all the requirements applied. 

b) The polar angle distribution of the Dt Dl~ combinations. 

c) The mass spectrum of the -D^-Ds~7isr combinations after 
subtraction of Ec.m. in case (2). d) The obtained M^+^.- 
spectrum. Cross-hached histograms show the normalized 
Mj^+ and Mj^t- sidebands contributions. The small con- 
tamination from the D*'^D*~ final state is shown by the open 
histograms. The signal windows are shown by vertical dashed 
lines. 

The Mf^^^.^{D*+D*-) distribution for M^.+^.- < 
5.0 GeV/c^ after all the requirements applied is shown in 
Fig.[3]a). A peak corresponding to e"'"e~ —5- D*"'"D*~7isr 
is again evident around zero recoil mass. We define 
the signal region for M^ecoiii^s~^ ) by a requirement 
±0.7( GeV/c^)^ around zero recoil mass. The polar an- 
gle distribution of the £'*"'"£'* ~ combinations and the 
mass spectrum of the -D*^-D*^7isr combinations (after 
subtraction of i?c.m.) that survive the M^^^^■^{D*+ D*~) 
requirement in case (2) are shown in Figs. [3]b,c). The 
full Mjyt+jjt~ spectrum after all the requirements are ap- 
plied is shown in Fig. [3]d). With such limited statistics 
no structure are evident. 

The contribution of multiple entries after all the re- 
quirements is found to be ~ 6%, ^ 22%, ^ 23%, for 
the D+D-, D+D*- and the D*+D*- final states, re- 
spectively. In such cases, only one D^*''^ D^*'' com- 
bination per event, that with the minimum value of 
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FIG. 3: a) The distributions of the M,2„^„ii(D*+D*" ) for 
M^, + ^,- < 5.0GeV/c^ after all the requirements applied. 

b) The polar angle distribution of the D*'^ D*~ combinations. 

c) The mass spectrum of the D^'^ D^~^isr combinations after 
subtraction of -Ec.m. in case (2). d) The M^.+^*- spectrum 
after all the requirements applied. Cross-hached histograms 



show the normalized M 



and M 



sidebands contribu- 
tions. The signal windows are shown by vertical dashed lines. 



used, where xl^j^+y xl,^j,-y xL(z3;+; 



cor- 



respond to the mass fits for the Dj, , and the 
D*~ candidates. 

The following sources of background arc con- 
sidered: (1) combinatorial background under the 
ni*'''^ (oi*'' ) peak combined with a correctly recon- 
structed ui*'* (oi*'''^) from the signal or other processes; 
(2) both the ui*-*^ and the Z?!*'' are combinatorial; (3) 
for the DfD^ final state: reflections from the e+e~ ^• 
D+D*^jisn and e+e~ £)*+£'*"7isr processes fol- 
lowed by D* — > I?s7, where the low momentum 7 is not 
reconstructed; for the D'^D*~ final state: reflection from 
the e+e~ D*+D*~7isr process followed by Z?* 
Ds^, where the low momentum 7 is not reconstructed; 
(4) reflection from the e+e" ^ L»^*^+Di*^"7r^,,,7isR 
processes, where the tt^Jj^^^ is not reconstructed. (5) the 
contribution of e+e^ — > D^s"^^ Di*'^~ ir^ , when an ener- 
getic tt" is misidcntificd as a single 7isr- 

The contribution from background (1) is extracted us- 
ing the ui*'" and D^^'^^ sidebands. Background (2) is 
present in both the M^(.)+ and M^(.)- sidebands and 
is, thus, subtracted twice. To account for this over- 
subtraction we use a double sideband region, when events 
are selected from both the M^(.)+ and the M^(.)- side- 
bands. The contributions from the combinatorial back- 
grounds (1-2) are shown in Figs. [T]d) and[2]d) as cross- 
hatched histograms. 

Backgrounds (3-4) arc suppressed by the tight re- 
quirement on A/j^j,^Qj[(Z?i*''^Z3i*''~). The remaining back- 
ground (3) for the D'^D*^ final state is estimated us- 
ing a MC simulation of the e+e^ — ?> D*~^D*~jisYi pro- 
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cess. To reproduce the shape of the D'^D*~ mass spec- 
trum we use the initial measurement of the D*~^D*~ 
mass spectrum. The remainder of background (3) for 
the DfD~ final state is estimated using a MC simulation 
of the e+e" D+D*^-fisn and e+e" L»*+L»*"7isr 
processes. To reproduce the shape of the DfD^ mass 
spectrum we use the initial measurement of the D*^ D*^ 
mass spectrum and the first iteration of the Df D*^ mass 
spectrum. The contributions from background (3) for 
the DfD^ and the D'^D*~ final states are shown in 
Figs. [T]a), [T]d) and Figs. [2]a), [2]d) as open histograms. 
Uncertainties in these estimates are included in the sys- 
tematic errors. 

To estimate the contribution from background (4), we 



7r°7iSR processes using 



study the 

fully reconstructed final states. From a MC study we es- 
timate the fraction of reconstructed events for the cases 
where the 7r° is not detected. After the application of the 
requirement on M'^^^^^^{D^s'^~^ D^s^ ) this contribution is 
found to be less than 0.5% and negligibly small; uncer- 
tainties in this estimate are included in the systematic 
errors. 

The contribution from background (5), in which an 
energetic Tr** is misidcntificd as the 7isr candidate, is de- 
termined from the data using fully reconstructed e+e~ — >■ 

< 



oi*^^ oi*^ tt" events. Only three events with M 



djd; 



the 7r° reconstruction efficiency. For the DfD*~ and the 
D*~^D*~ final states the expected backgrounds are ^ 0.6 
events and events in the whole Mj^+jj,- and M^.+^.- 
mass ranges. The probability of tt'' 7 misidentification 
due to asymmetric tt*^ — > 77 decays is also estimated to 
be small. Thus the contribution from background (5) is 
found to be negligibly small; uncertainties in these esti- 
mates are included in the systematic error. 

The e+e~ — > Di*''^ Di*''~ cross sections are extracted 
from the background subtracted Z?i*''^Z?i*'' mass distri- 
butions 



dN/dm 
ritotdL/dm ' 



(2) 



where m = M 



dN/dm is the obtained mass 
spectrum, rjtot is the total efficiency and the factor 
dL/dm is the differential ISR luminosity [l^. The to- 
tal efficiencies determined by the MC simulation grow 
quadratically with energy from O.Of 5%, 0.010%, 0.005% 
near threshold to 0.045%, 0.025%, 0.011% at 5.0 GeV/c^ 
for the D+DJ, D+ D*' and the D*+ D*' final states, 

respectively. The resulting e+e^ Dl Dl exclu- 
sive cross sections averaged over the bin width are shown 
in Fig. 21 Since the bin width is much larger than the 
7\/ (.)+ (.)_ resolution, which varies from ^ 2MeV/c^ 



5.0GeV/c2 and M^+^-^o - ^c.m 
in the data. Assuming a uniform tt" polar angle distribu- 
tion, this background contribution in the | cos(0^+^- ) | > 
0.9 signal sub-sample (case 1) is 3events/9e^o ~ 0.6 
events in the whole M^+^- mass range, where t-^o is 



> 0.5 GeV are found around threshold to 



6MeV/c2 at M 



5.0GcV/c^, no correction for resolution is applied. The 
next-to-leading order radiative corrections are taken into 
account by the dL/dm formula. The next-to-next-to- 
leading order corrections are included in the systematics. 
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FIG. 4: The cross section averaged over the bin width for a) the e^e -^DfD^ process; b) the e"'"e -^DfD* + cc. process; 
c) the e+e" D'+ Dl~ process. Error bars show statistical uncertainties only. There is a common systematic uncertainty 
for all measurements, 11% for DfD~, 17% for DfD*~ and 31% for D*^D*~ . This uncertainty is described in the text The 
dotted hnes show masses of the V(4040), VCllSO) and V(4415) states [l3 |. 



The R ratio, defined as i? = (j{e^e^ — > 
hadrons) / a{e^ — > where cr(e^e~ fi^fi^) = 

Atto^/Ss, for the sum of the exclusive e^e~ 

cross sections is shown in Fig. \b\ 

The systematic errors for the a{e'^e~ — ?► Z?i*''^Z?i*'' ) 
measurements are summarized in Table ID The system- 
atic errors associated with the background (1-2) subtrac- 
tion are estimated from the uncertainty in the scaling fac- 



tors for the sideband subtractions. This is done using fits 
to the M^(«)+ and M^(,)- distributions in the data with 
different signal and background parameterizations and 
are found to be 3%, 7% and 24% for the D+ , D+D*- 
and the Dl^Dl" final states, respectively. Uncertainties 
in the contribution from background (3) are estimated 
to be 2% for the DfDj final state and smaller than 1% 
for the DfD*~ final state. Uncertainties in the back- 
grounds (4-5) are estimated conservatively to be smaller 



6 



0.3 
0.2 
0.1 




-\ 



X4 



h 



++ 



3.8 



4.2 4.4 4.6 4.8 f 

M(D (*)+D (*)-) GeV/c^ 



FIG. 5: The R ratio for the sum of the Df D~ , Df D*~ and 
the D*'^ D*~ final states. The vertical dotted lines show the 
masses of the ^'(4040), i/)(4160) and V(4415) states [H]. Error 
bars show statistical uncertainties only. 



TABLE I: Contributions to the systematic error on the cross 
sections. 



Source 


DtD- 


Dtor 


d*+d:- 


Background subtraction 


±4% 


±7% 


±24% 


Cross section calculation 


±7% 


±ff% 


±12% 




±5% 


±5% 


±5% 


Reconstruction 


±6% 


±fO% 


±16% 


Kaon identification 


±2% 


±2% 


±2% 


Total 


±ff% 


±f7% 


±31% 



than 1% of the signal for all final states. The systematic 
errors ascribed to the cross section calculation include an 
error on the differential ISR luminosity (2%) and statis- 
tical errors of the MC simulation for the total efficiency 
calculations. In the case of the DfD*~ state there is 
an additional uncertainty due to the unknown helicity 
distribution, estimated following a procedure similar to 
that used in Ref Q . Another source of systematic error 
comes from the uncertainties in track and photon recon- 
struction efficiencies (1% per track, 1.5% per photon and 
7% per soft photon). Other contributions include the un- 
certainty in the identification efficiency and the absolute 
Df and D*~^ branching fractions [l^l . The total system- 
atic uncertainties are 11%, 17% and 31% for the DfD'^, 
D+D*~ and the Dl+D*~ final states, respectively. 

In summary, we report the measurement of the 
e+e- DV'^DV' exclusive cross sections over the 
center-of-mass energy range from the Ds thresh- 
olds to 5.0 GeV. A clear peak at threshold around the 
■0(4040) mass is seen in the e'*"e~ — ^ DfD~ cross sec- 
tion. In the e+e~ D^D*~ cross section two peaks 
are evident around the '0(4160) and the -0(4415) masses. 



The limited statistics do not reveal any structures in the 
e+e^ — !• D*^Dl~ cross section. The obtained R ratio for 

the sum of e^e~ Dl Dl cross sections has a rich 
structure including peaks around the ^^(4040), -0(4160) 
and the 0(4415) masses. Both the e^e" — > DfD*~ cross 
section and the R ratio exhibit an obvious dip near the 
y(4260) mass, similar to what is seen in e'*"e~ — > D*D* 
and in the total cross section for charm production. The 
obtained cross sections are consistent within errors with 
those from BaBar [T^]. The CLEO exclusive cross sec- 
tions are not directly comparable to those from Belle 
as they are not radiatively corrected, but generally seem 
to reflect consistency. 

In this study we do not perform a fit to the obtained 
e+e^ — > cross sections. The numerous open 

charm thresholds in this region complicate the cross sec- 
tions behavior and coupled-channel modifications to the 
description of any particular final state require one to 
take into account all other final states contributing to 
the total cross section for charm production. 
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Recently a number of measurements of exclusive cross 
sections for e'^e~ annihilation into charmed hadron pairs 
above open-charm threshold were performed by the B- 
factory experiments using initial-state radiation (ISR). 
These include Belle measurements [l| of e+e" cross sec- 
tions for the DD {D = D° or D+), D+D*-, D*+D*-, 
D°D-Tr+, L>°£)*-7r+ and A+A7_final_states [|Ji and 
BaBar measurements of the DD, DD*, D*D* final 
states 0,11], which arc, in general, consistent with those 
of Belle. In addition, CLEO scanned the e~^e~ center 
of mass energy range from 3.97 to 4.26 GeV and ob- 
tained exclusive cross sections for the DD, DD* , D*D* , 
DDtt and DD*tt final states @. The first measure- 
ments of the exclusive cross sections for e+e" annihila- 
tion into charmed strange meson pairs Ds Ds were 
performed by CLEO with high accuracy but with limited 
maximum energy (4.26 GeV) ^. Recently BaBar pre- 

sented e+e" Dl '^Dl ' cross sections averaged over 
100 MeV wide bins [l^l ■ The observed cross sections were 
found to be an order of magnitude smaller than those for 
non-strange charmed meson production. 

Although the recent BES fit to the total cross section 
for hadron production in e+e~ provided new parame- 
ter values for the ip resonances [11|, the available exclu- 
sive e'^e~ cross sections have not yet been qualitatively 
explained. One of the main problems is the numerous 
open charm thresholds in the region that influence the 
cross section behavior and, thus, complicate theoretical 
descriptions. 

The Y states [l^ (with masses above open charm 



threshold and quantum numbers J^'~' = 1 ), which do 
not exhibit strong decays to any of the measured open 
charm final states and have remain unexplained since 
their discovery more than five years, provide additional 
motivation to pursue all possible experimental informa- 
tion about the decomposition of charmed particle pro- 
duction in the charm-threshold region. 

Here we report a measurement of exclusive e+e~ — >^ 
£'i*'''''Z?i*'' cross sections as a function of ccnter-of-mass 
energy from the D), ' Dl ' thresholds to 5.0 GeV, con- 
tinuing our studies of the exclusive open charm produc- 
tion in this mass range. The analysis is based on a study 
of events with ISR photons in a data sample collected 
with the Belle detector [l^ at the T{nS) (n = 1,...,5) 
resonances and nearby continuum with an integrated 
luminosity of 967 fb"^ at the KERB [13 asymmetric- 
energy e'^e" collider. 

We follow the full reconstruction method that was pre- 
viously used for the measurements of the e'^e' cross sec- 
tions to D'D, D°D-TT+ and D°D*-'k+ final states 0, 

13, [1]. We select e+e^ D^*^^ Di*^^jisji signal events 
in which the Ds and Ds ' mesons are fully recon- 
structed. ISR photon candidates are indicated by 7isr. 
In general, an 7isr is not required to be detected and its 
presence in the event is inferred from a peak at zero in the 
spectrum of recoil mass squared against the Ds Ds 
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system. The recoil mass squared is defined as: 

ACo.,(i^(*)+i?W-) = - _ 

where Ec.m. is the initial e+e^ ccnter-of-mass (cm.) 
energy, E^(,)+^m- and are energy and 

three- momentum of the Z^s Dg combination, respec- 
tively. To suppress backgrounds two cases are consid- 
ered: (1) the 71SR is outside of the detector acceptance 
and the polar angle for the Ds Ds combination 
in the cm. frame is in the range | cos(0^(.)+^(.)-)| > 
0.9; (2) the 71SR is within the detector acceptance 
(I cos(0^(.)+^(.)- )| < 0.9). In the latter case, the 
71SR is required to be detected and the mass of the 
71SR combination should to be greater than 
{Ec.m. — 0.5 GeV). To suppress backgrounds from 

e+e- ^ Di*)+Di*'-(m)(7r+7r-)7isR,(m = 1,2,...) pro- 
cesses, we exclude events that contain additional charged 
tracks that were not used in the Ds and the Dg re- 
construction. 

All charged tracks are required to originate from the 
vicinity of the interaction point (IP); we impose the re- 
quirements \dr\ < 1cm and \dz\ < 4cm, where dr and 
dz are the impact parameters perpendicular to and along 
the beam direction with respect to the IP, respectively. 
Charged kaons are required to have a ratio of particle 
identification likelihood, Vk = jO-k /{jO-k + jO-tt) > 0.6 (l5| . 
No identification requirements arc applied for pion can- 
didates. 

Kg candidates are reconstructed from tt+tt" pairs with 
an invariant mass within 10 MeV/c^ of the Kg mass. The 
distance between the two pion tracks at the Kg vertex 
must be less than 1 cm, the transverse flight distance 
from the interaction point is required to be greater than 
0.1 cm, and the angle between the Kg momentum direc- 
tion and the flight direction in the x — y plane should be 
smaller than 0.1 rad. 

Photons arc reconstructed in the electromagnetic 
calorimeter as showers with energies greater than 50 MeV 
that are not associated with charged tracks. Pairs of pho- 
tons are combined to form 7r° candidates. If the mass of 
a 77 pair lies within 15MeV/c^ of the tt" mass, the pair 
is fltted with a tt^ mass constraint and considered as a tt" 
candidate. ISR photon candidates are required to have 
energies greater than 2.5 GeV. Photon candidates used 
in ?7, rj' and D*~^ reconstruction are required to have en- 
ergies greater than 100 MeV. 

rj candidates are reconstructed using Tr+Tr^Tr" 
(±10MeV/c2 mass window) and 77 (±20MeV/c2 
mass window) decay modes (~ 2.5 cr in each case), rj' 
candidates are reconstructed using t^tt+tt" (±10MeV/c^ 
mass window) and jtt^tt^ (±15McV/c^ mass window) 
decay modes (~ 2.0 u in each case). A mass- and 
vertex-constrained fit is applied to r/ and rj' candidates. 

candidates are reconstructed using six decay 
modes: + , K-K+tt+, K-K+tt+tt^, KIK-t:+tt+, 



T]TT^ and t]'tt^. Before calculation of the D'^ candidate 
mass, a vertex fit to a common vertex is performed for 
tracks that form the Df candidate. A ±15 MeV/c^ mass 
signal window is used for all modes (^ 3ct in each case). 
To improve the momentum resolution of Df meson can- 
didates, the tracks from the Df candidate are fitted to 
a common vertex with a mass Df mass constraint. D*'^ 
candidates are reconstructed using the Dfj decay mode. 
A ±15MeV/c^ mass window is used (^ 2.5 a). A mass- 
constrained fit is also applied to D*"*" candidates. 

The Df and D*~^ sidebands used for background stud- 
ies are four times as large as the signal region and are di- 
vided into windows of the same width as that of the sig- 
nal. To avoid signal over-subtraction, the selected side- 
bands are shifted by 30MeV/c^ from the signal region. 
The Ds{D*) candidates from these sidebands are refitted 
to the central mass value of each window. 

The Mrocoii(^.t^r) distribution for Mj^^^- < 
5.0GeV/c^ after all the requirements described above 
is shown in Fig. [T]a). A clear peak corresponding to 
the e+e" D+-D771SR process is evident near zero 
recoil mass. The shoulder at positive values is mainly 
due to e+e~ — > DfD*~'yisn background. We de- 
fine the signal region with an asymmetric requirement 
-0.7(GcV/c2)2 < Ml^^,,{DtDj) < 0.4(GeV/c2)2 to 
suppress event this background. The polar angle distri- 
bution of the DfDj^ combinations and the mass spec- 
trum of the Df DjjisK combinations (after subtraction 
of Ec.m.) in case (2), after M^^^^^^{Df D^) requirement, 
arc shown in Figs. [T]b,c). These distributions are in 
agreement with a MC simulation and are typical of ISR 
production. The Af^+^- spectrum with all the require- 
ments applied is shown in Fig.[l]d). A clear peak is seen 
at threshold near the -0(4040) mass. 

The M^^^^ii (£)+£)*-) distribution for A/^+^.- < 
5.0 GeV/c^ after all the requirements arc applied is shown 
in Fig. [2]a). A clear peak corresponding to e+e" 
£'+D*~7isR signal process is again evident around zero. 
We define the signal region for -^recoil (-^s"-^*") by a re- 
quirement ±0.7( GeV/c^)^ around zero. The polar angle 
distribution of the DfD*^ combinations and the mass 
spectrum of D+D*~7isr combinations (after subtraction 
of the -Ec.m. energy) in case (2) after the requirement on 
^•^recoii(-D^£T) is apphed are shown in Figs.[l]b, c). The 
Mjj+jjt- spectrum after all the requirements are applied 
is shown in Fig. [2]d). Two clear peaks are seen at the 
■0(4160) and the '!/'(4415) masses. 

The A/4^^ii(£)*+£);-) distribution for < 
5.0 GeV/c^ after all the requirements applied is shown in 
Fig.[3]a). A peak corresponding to e+e^ — > D*^D*~ ^isji 
is again evident around zero recoil mass. We define 
the signal region for M^ccoiii^t^ ^t ) by a requirement 
±0.7( GeV/c^)^ around zero recoil mass. The polar an- 
gle distribution of the D*~^D*~ combinations and the 
mass spectrum of the D*"'"D*~7isr combinations (after 
subtraction of Ec.m.) that survive the M^^^^-^{D*+ Dl~) 
requirement in case (2) are shown in Figs. [3]b, c). The 
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FIG. 1: a) The distribution oi Mf^^^^iiDj D~) for M^+^- < 

5.0GeV/c^ after all the requirements are applied, b) The 
polar angle distribution of the DfD~ combinations, c) The 
mass spectrum of the Df D^jisn combinations after subtrac- 
tion of Ec.m. energy in case (2). d) The M^+^- spectrum 
after all the requirements applied. Cross-hached histograms 
show the normalized and M^- sideband contributions. 

Feed down from the D^D*^ final state is shown by the open 
histograms. The signal windows are shown by vertical dashed 
lines. 



FIG. 2: a) The distribution of the M^^^^n{Dt D*~ ) for 
Mj^+j^,- < 5.0GeV/c^ after all the requirements applied. 

b) The polar angle distribution of the DfD*~ combinations. 

c) The mass spectrum of the DfO^'jisii combinations after 
subtraction of Ec.m. in case (2). d) The obtained M^+^.- 
spectrum. Cross-hached histograms show the normalized 

and M^t^ sidebands contributions. The small con- 
tamination from the D^'^Dg^ final state is shown by the open 
histograms. The signal windows are shown by vertical dashed 
lines. 



full Mjj,+jj,- spectrum after all the requirements are ap- 
plied is shown in Fig. [3]d). With such limited statistics 
no structures are evident. 

The contribution of multiple entries after all the re- 
quirements is found to be ^ 6%, ^ 22%, ~ 23%, for 
the D+Dj, D+D*- and the D*+-D*~ final states, re- 
spectively. In such cases, only one D^*''^ D^*'' com- 
bination per event, that with the minimum value of 

Xtot — ^M{Dt) ~^ ^M{D7) ~^ ^^M{Dt + )'^ ~^ ^^M(Dl-)^'' '^^ 

cor- 



used, where xlj^j^+y xl,^j,~y xL(d*+) ^""^ AiiDl , 
respond to the mass fits for the _D*+ and the 

D*~ candidates. 

The following sources of background are con- 
sidered: (1) combinatorial background under the 



peak combined with a correctly recon- 



structed Z?i*-'~(_Di*-'^) from the signal or other processes; 

(2) both the Di*"*^ and the D^*'' are combinatorial; (3) 
for the DfD^ final state: reflections from the e+e^ — s- 
DfD*~'-fisR and e'^e~ D*+D*~7isr processes fol- 
lowed by D* — > £'s7, where the low momentum 7 is not 
reconstructed; for the D^D*~ final state: reflection from 
the e+e~ — >• £'*+£'*~7isr process followed by D* 
Dsj, where the low momentum 7 is not reconstructed; 



(4) reflection from the 



-> D. 



processes, where the ir^iss ^'^^ reconstructed. (5) the 



contribution of e^e 



D 



(*)+ n(*) 



Di*>~Tr'^, when an ener- 



getic TT^ is misidentified as a single 7isr. 

The contribution from background (1) is extracted us- 



ing the oi*^ and D^*'^ sidebands. Background (2) is 
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3: a) The distributions of the M^^coiiiDs^ D*' ) for 
^* + ^,_ < 5.0GeV/c^ after all the requirements applied. 

b) The polar angle distribution of the D^'^ Dt," combinations. 

c) The mass spectrum of the Dg'^D^~'yisn. combinations after 



subtraction of -Ec.m. in case (2). d) The Mj- 



spectrum 



after all the requirements applied. Cross-hached histograms 



show the normalized M 



and M 



sidebands contribu- 
tions. The signal windows are shown by vertical dashed lines. 



71SR present in both the M 



and M, 



sidebands and 



is, thus, subtracted twice. To account for this over- 
subtraction we use a double sideband region, when events 
are selected from both the M^(.)+ and the M^(.)- side- 
bands. The contributions from the combinatorial back- 
grounds (1-2) are shown in Figs. [T]d) and[2]d) as cross- 
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hatched histograms. 

Backgrounds (3-4) are suppressed by the tight re- 
quirement on M^^j,j.Q;|(£>i*-'^£'i*^ ). The remaining back- 
ground (3) for the D'^D*~ final state is estimated us- 
ing a MC simulation of the e+e^ -D*+L'*~7isr pro- 
cess. To reproduce the shape of the DfDl^ mass spec- 
trum we use the initial measurement of the D*^D*~ 
mass spectrum. The remainder of background (3) for 
the DfD~ final state is estimated using a MC simulation 
of the e+e~ £>+I?*~7isr and e+e~ D*+Dl~ jisn 
processes. To reproduce the shape of the DfD~ mass 
spectrum wc use the initial measurement of the D*^D*~ 
mass spectrum and the first iteration of the D'^D*~ mass 
spectrum. The contributions from background (3) for 
the DfDj^ and the D'^D*~ final states are shown in 
Figs. [T]a), [T]d) and Figs. [Ha), [5]d) as open histograms. 
Uncertainties in these estimates are included in the sys- 
tematic errors. 

To estimate the contribution from background (4), we 

study the e+e^ — s> Di*^~^ Di*^~ tt^jisr processes using 
fully reconstructed final states. From a MC study we es- 
timate the fraction of reconstructed events for the cases 
where the 7r° is not detected. After the application of the 
requirement on M^^^^-^-i{Di*'^~^ oi*"^ ) this contribution is 
found to be less than 0.5% and negligibly small; uncer- 
tainties in this estimate are included in the systematic 
errors. 

The contribution from background (5), in which an 
energetic tt" is misidentified as the 7isr candidate, is de- 
termined from the data using fully reconstructed e+e~ 
Di*^^Di*^ tt" events. Only three events with Mj^+j^- < 
5.0GeV/c^ and Mj^+j^-^o - -Bc.m. > 0.5 GeV are found 
in the data. Assuming a uniform tt" polar angle distribu- 
tion, this background contribution in the | cos(0^+^- ) | > 
0.9 signal sub-sample (case 1) is 3events/9e^o ~ 0.6 
events in the whole M 



DtD: 



mass range, where €^^0 is 
the tt" reconstruction efficiency. For the DfD*~ and the 
final states the expected backgrounds are ~ 0.6 



events and events in the whole Mj^+j^-,- and Mj^^+j^*- 
mass ranges. The probability of tt" — > 7 misidentification 
due to asymmetric 7r° — >■ 77 decays is also estimated to 
be small. Thus the contribution from background (5) is 
found to be negligibly small; uncertainties in these esti- 
mates are included in the systematic error. 

cross sections are extracted 
mass distri- 



from the background subtracted Z) 
butions 



a{e+e 



where m = M 



dN/dm 
rjtotdL/dm ' 



(2) 



dN/dm is the obtained mass 



spectrum, ?7tot is the total efficiency [17 
dL/dm is the differential ISR luminosity: 



The factor 



dL/dm = — ( (2 - 2a: 
where x — I 



I In 



1 + C 
l-C 



2mC 



,(3) 



m^/£'c.m.^, C is the total integrated lu- 
minosity and C = cos 00, where 0o denotes the polar 
angle range for 71SR in the e^e~ cm. frame: < 
^7isR < ~ ^0- The total efficiencies determined 
by the MC simulation grow quadratically with energy 
from 0.015%, 0.010%, 0.005% near threshold to 0.045%, 
0.025%, 0.011% at 5.0GeV/c2 for the D+ Dj , D+D*- 
and the D*^ D*~ final states, respectively. The result- 
ing e+e~ 

exclusive cross sections aver- 
aged over the bin width are shown in Fig. 3] Since the 
bin width is much larger than the M^(.)+^(.)- resolu- 
tion, which varies from ~ 2MeV/c^ around threshold to 
6MeV/c2 at M^(.)+^(.)- = 5.0GeV/c2, no correction 
for resolution is applied. The ncxt-to-leading order ra- 
diative corrections are taken into account by the dL/dm 
formula. The next-to-next-to-leading order corrections 
are included in the systematics. The contribution of fi- 
nal state radiation (FSR) is strongly suppressed [l^ and 
is neglected in this study. 



gO.6 



0.4 
0.2 



„Lipv4 



a) 



3.8 4 4.2 4.4 4.6 4.8 

M(D +D -), GeV/c^ 




4.2 4.4 4.6 4.8 5 

M(D +D *-), Ge\/c 



M(D *+D *-), Ge\/c 



FIG. 4: The cross section averaged over the bin width for a) the e^e -^DfD^ process; b) the e'''e -^DfD* -I- c.c. process; 
c) the e+e" ^ D*+D*' process. Error bars show statistical uncertainties only. There is a common systematic uncertainty 
for all measurements, 11% for Dj , 17% for DfD*~ and 31% for D*^D*~ . This uncertainty is described in the text. The 
dotted hnes show masses of the V(4040), 1^(4160) and 7/)(4415) states [la |. 



The R ratio, defined as R 



(^(e+e 



hadrons) / a{e^ e — > /i^/i ), where a[e^e /^^M ) 
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Ana'^/Ss, for the sum of the exclusive e+e 
cross sections is shown in Fig. [5j 




4.2 4.4 4.6 

M(D (*)+D (*)-) 



FIG. 5: The R ratio for tlie sum of the Df , Df D*' and 
the 0*^0*" final states. The vertical dotted lines show the 
masses of the «/;(4040), -0(4160) and V'(4415) states [IB]. Error 
bars show statistical uncertainties only. 

The systematic errors for the (T(e+e^ oi*^^ d[*'^ ) 
measurements are summarized in Table HI The system- 



TABLE I: Contributions to the systematic error on the cross 
sections. 



Source 


DtD- 


DtDT 


D*+Dr 


Background subtraction 


±4% 


±7% 


±24% 


Cross section calculation 


±7% 


±11% 


±12% 




±5% 


±5% 


±5% 


Reconstruction 


±6% 


±10% 


±16% 


Kaon identification 


±2% 


±2% 


±2% 


Total 


±11% 


±17% 


±31% 



atic errors associated with the background (1-2) subtrac- 
tion are estimated from the uncertainty in the scaling fac- 
tors for the sideband subtractions. This is done using fits 
to the M^(*)+ and M^(,)- distributions in the data with 
different signal and background parameterizations and 
are found to be 3%, 7% and 24% for the D+ D- , D+D*- 
and the D*^D*~ final states, respectively. Uncertainties 
in the contribution from background (3) are estimated 
to be 2% for the DfD~ final state and smaller than 1% 
for the DfD*~ final state. Uncertainties in the back- 
grounds (4-5) are estimated conservatively to be smaller 
than 1% of the signal for all final states. The systematic 
errors ascribed to the cross section calculation include an 
error on the differential ISR luminosity (2%) and statis- 
tical errors of the MC simulation for the total efficiency 
calculations. In the case of the D+_D*~ state there is 
an additional uncertainty due to the unknown helicity 
distribution, estimated following a procedure similar to 
that used in Ref Q. Another source of systematic error 
comes from the uncertainties in track and photon recon- 
struction efficiencies (1% per track, 1.5% per photon and 
7% per soft photon). Other contributions include the un- 
certainty in the identification efficiency and the absolute 



Df and D*~^ branching fractions [la| . The total system- 
atic uncertainties are 11%, 17% and 31% for the DfD~, 
D+D*~ and the D*+D*^ final states, respectively. 

In summary, we report the measurement of the 
— >■ ni*'''^ oi*^^ exclusive cross sections over the 



center-of-mass energy range from the Dl '^Dl ' thresh- 
olds to 5.0 GeV. A clear peak at threshold around the 
-0(4040) mass is seen in the e+e^ — ?> DfDj cross sec- 
tion. In the e+e" — t- DfDl~ cross section two peaks 
are evident around the -0(4160) and the ?/;(4415) masses. 
The limited statistics do not reveal any structures in the 
e+e^ Dl~^D*~ cross section. The obtained R ratio for 
the sum of e+e^ — > D), cross sections has a rich 

structure including peaks around the -0(4040), 0(4160) 
and the -0(4415) masses. Both the e+e^ — > DfDl~ cross 
section and the R ratio exhibit an obvious dip near the 
F(4260) mass, similar to what is seen in e~^e~ — > D*D* 
and in the total cross section for charm production. The 
obtained cross sections are consistent within errors with 
those from BaBar [l^. The CLEO exclusive cross sec- 
tions are not radiatively corrected and, therefore, can- 
not be directly compared to the results reported here. 

In this study we do not perform a fit to the obtained 
e+e^ — > Ds Ds cross sections. The numerous open 
charm thresholds in this region complicate the cross sec- 
tions behavior and coupled-channel modifications to the 
description of any particular final state require one to 
take into account all other final states contributing to 
the total cross section for charm production. 
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